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Abstract 

It is argued that the dynamics of an isolated system, due to the concrete procedure by 
which it is separated from the environment, has a non-Hamiltonian contribution. By a 
unified quantum field theoretical treatment of typical subdynamics, e.g., hydrodynamics, 
kinetic theory, master equation for a particle interacting with matter, we look for the 
structure of this more general dynamics. 

> 

O : 1 The concept of physical system 

q \ Quantum mechanics (QM) has non-separability as its most striking feature, i.e., one can- 
not attribute "properties" to parts of a system and therefore typical problems like the 
measurement process and EPR situations arise. This feature is so deeply rooted in the 
mathematical structure of QM that we believe one should not try to make it less strin- 
gent, e.g., by attempts like "spontaneous reduction" We prefer instead to weaken the 
very concept of physical system: usually the "isolation" of a physical system is taken for 
granted, while in our opinion the way in which isolation is achieved belongs to the very 
definition of the system. Any attempt inside QM to obtain the subdynamics for a subsys- 
tem enforces the introduction of a suitable time scale in order to break the correlations 
with the environment; in a completely sharp description of the dynamics of a subsystem 
the physics of the whole universe would enter. The preparation procedure leading to a 
system, isolated during a time interval \to,t\\ and confined in a spatial region u, covers a 
time interval [T, to], that will be called "preparation time". Due to the confinement the 
basic space-time symmetries are broken and by suitable boundary conditions "peculiar" 
properties of the system are introduced. This obviously reduces the universal character of 
the dynamical description, however an important universal behaviour still remains due to 
symmetry and locality (or short range character) of effective interactions, whose relevance 
becomes particularly evident in the quantum field theoretical approach. We thus regard 
a physical system as a part of the world under control by a suitable preparation, whose 
local behaviour is explained in terms of locally interacting quantum fields. The choice of 
these fields depends on the level of description of the system. A large part of physics can 
be explained in terms of quantum fields related to molecules with a typical time scale of 
the order pa 10~ 13 s; a much more refined description arises if the basic fields are related to 
nuclei and electrons, then the basic theory would be QED and a much smaller time scale 
pa 10~ 23 s could be considered: however such role of QED as basic theory of macrosystems 
is far from being exploited. 

In a sense our viewpoint can appear as opposite to the most widely spread one, which we 
synthetise as follows: particles are the primary systems, related to non-confined quantum 
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fields and to basic symmetries, all other systems are structure of particles; one then tries 
to obtain a typical macroscopic behaviour in some suitable thermodynamic limit. Accord- 
ing to us, on the contrary, macroscopic systems are to be taken as the primary systems, 
even if in their definition time scales and spatial confinement must be carefully taken into 
account; the theoretical framework for their description is quantum field theory, locality 
and quantisation taking the place of the atomistic model. In this context particles are 
a derived concept. The description of non-equilibrium systems is put in the foreground 
and at least in principle should be performed taking boundary effects into account; pro- 
cedures like "continuous" limit should be applied only at the end, if one wants to get rid 
of boundary effects. This standpoint is closer to thermodynamics and electromagnetism, 
while the former one originates from classical mechanics. The relevance of macroscopic 
systems for the foundations of QM is the starting point of Ludwig's axiomatic approach 
to QM [g. The insistence on the distinction between these two attitudes is due to the 
fact that they lead in a natural way to two different formulations of the dynamics. In 
the first approach one associates a wave function ip to each system (^(x, t) for one par- 
ticle, . . . , ■0(xi,x 2 , . . . , Xjv,£) for N particles); obviously if one describes situations like 
"unpolarised" particles it is appropriate to use a statistical operator, in order to take 
a lack of control of the experimental specification into account. This aspect becomes 
increasingly important for large N, so that statistical operators are very useful for macro- 
scopic systems, nevertheless the basic dynamics is given by an evolution operator for the 
wavefunction ip. On the other hand, starting with a macroscopic system, one is led to 
assume a statistical operator g t as the most appropriate mathematical representation of 
the preparation procedure until time t. The set K,(7i) of statistical operators on the 
Hilbert space TL becomes most important and the space T{TL) of trace-class operators, 
which is generated in a natural way by ICiTi) [fC(TL) is the base of the base-norm space 
T(7~0], plays a role similar to that of 7i in the previous formalism. Correspondingly 
unitary operators on 7i in the first approach are replaced by affine maps of K,(7i) in 
/C(7i), i.e. by positive, trace-preserving maps on T{7i). If the system is isolated in the 
time interval [to^i] the spontaneous repreparations g t , t G [to,ti], are related together by 
Qt = M-tt'Qt' (t > t'), where the evolution system {M. t "t' t" > f} satisfies the composition 
rule M.t"'t' — M-t'"t"M-t"t' (t' — t" < t'"). We stress the fact that there is no reason 
to assume that /4ft' has an inverse. If M.J,, t , exists then /Aft' = Uft' • Uft> ( see f° r 
example ||), with lift' unitary or antiunitary operator and one is brought back to the 
Hilbert space formalism: tpf = Uft'^Pt 1 - The dynamics in the present framework is indeed 
more general and has irreversibility as typical phenomenon. 

To determine the maps M.t"t' a choice of relevant observables is necessary: when a 
time scale is introduced, only those observables should be considered, whose expectation 
values do not appreciably vary in a time interval of the order of the time scale. It is 
thus necessary to work in the Heisenberg picture, i.e. with the adjoint map /A tto , and 
consider expressions of the form Ai' tto A, A being a relevant observable. For the same 
system different descriptions can be given by different choices of relevant observables 
and corresponding time scales: e.g., hydrodynamic or kinetic description of a continuum. 
Skipping questions of mathematical rigour we can assume the differential equation 

^M' tt0 = C' t M' tt0 , (1.1) 
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and represent Ai' tto in the form Ai' tto = T ye^o *'J in terms of the generator C' t . It 

is well known (rigorously for bounded C' t ) that if M.' t u t i has the additional property of 
complete positivity (CP), C t has the Lindblad structure Q: 

C' t B = + l - (H t B - BH t ) - ~ (A t B + BA t ) + l - £ L f tj BL t j (1.2) 

i 

In our framework the assumption of CP can appear too restrictive since only a suitable 
subset of observables is relevant and one expects that a modified concept of CP of M. t n t , 
relatively to these observables should be given, leading to a more general structure of C t ; 
we shall return to this point in the sequel. 

The more general description of the dynamics that we are considering allows the intro- 
duction of the concept of trajectory in quantum theory. In fact in the general formalism of 
continuous measurement approach || one has that an evolution system {A4t"t' t" — t'} 
with Ct having the Lindblad structure can be decomposed on a space Y^ of trajectories 
for stochastic variables £(£). The concept of trajectory is particularly useful in the field 
of quantum optics, as is shown by the many interesting applications to be found in the 
literature. More precisely one can consider Wiener and jump processes related to the 
operators L t j in the sense that the expectations of the increments for example in the case 
of jump processes are given by 0: 

(d^(t)) = dtTr((Ll J L tj )g(t))- 

One can define cx-algebras B(Y t t l ") of subsets Q\, of Y t f and construct operation valued 
measures T\, (fi*/ ) on B(Y^, ) in such a way that M. t "t' = ^(Xt 1 )• Then for any decom- 
position Y/j = U Q (f2*/ a ) with disjoint subsets VL\, a one has 

aw = £^:'j. (i.3) 

a 

One can therefore claim that the quantum dynamics of the system is compatible with 
the evolution of classical stochastic variables; typically the probability that the trajec- 
tory of these variables for t' < t < t" belongs to a subset f2|, a is given by p{Sl\> a ) = 
Tr (Q^,' a ) g(t')^j . The decomposition of A4t"f by operation valued stochastic processes 

T\, on a suitable trajectory space Yf is not unique, i.e. there are many compatible objec- 
tive "classical" pictures which are consistent with the quantum evolution, a feature that 
can be linked with a "generalised concept of complementarity". The very possibility of 
recovering some kind of classical insight into QM is due to the non-Hamiltonian evolution; 
obviously ( |1.3| ) would be inconsistent with M. t "t' = U t "t' • U}„ t ,, since for g t > = |t/V)(?/v| 
the l.h.s. of (|1.3j ) is a pure state and the r.h.s. is a mixture. 

In the framework we have now presented dynamics is given by C t 0; while one expects 
that the Hamiltonian part is fixed by the local interactions, the remaining part is con- 
nected to the preparation procedure of the isolated system: it cannot be strictly derived 
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from a Hamiltonian theory of a larger system surrounding it (in fact also for this larger 
system an C t should be determined). The problem arises to give reasonable criteria for the 
construction of C t . One can expect that near to equilibrium only the Hamiltonian part 
j?[H, •] is important, as it is clearly indicated by the great success of equilibrium statistical 
mechanics; however the non-Hamiltonian part of C t is relevant for irreversible behaviour 
and for a full explanation of approach to equilibrium: in fact energy conservation £ t H = 
(at least on the relevant time scale) grants for the existence of an "eigenstate" of C t with 
practically zero eigenvalue, while one expects that the other eigenvalues of C t have a 
negative real part. 



2 A microsystem interacting with matter 

We shall take on later the problem of an explicit construction of g(t) for a macrosystem 
Sm- We assume now that g(t) is given (e.g., the system is at equilibrium) and consider 
the problem of describing the new system S = Sm + a microsystem, S still being confined 
inside a region ou. The Hamiltonian of S is: 



H = H + H M + V # = E E f^f 

f 



a f , a £ 
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where Ef are the eigenvalues of the operator — |^A 2 : 

h 2 

~^ A 2^/(x) = E f Uf{x) u/(x) = xedu. 
Let us assume for the statistical operator of S the following structure: 



6 (t) = J2alQM(t)a f Q gf (t), 
gf 



(2.4) 



ajQu = 0. 



(2.5) 



The coefficients g g f build a positive, trace one matrix, which can be considered as the 
representative of a statistical operator g^(t) in the Hilbert space HS l > spanned by the 
states Uf (g g f = (uglg^'luf)). Equation ( |2.5|) indicates that the system Sm has charge 
Q = S/O/Q/ with value zero, i.e. it does not contain the microsystem. Equation ( |2.4| ) 
represents the fact that Sm has been perturbed by the additional particle and there- 
fore presents a new dynamical behaviour contained in the coefficients g g f{t) that can 
be picked out studying the time evolution of the observables A = J2hk ^h^hk^k] m fact 
the subdynamics of these observables provides the QM of a one-particle system with 

Hilbert space statistical operator g^ and observables A*' , with matrix elements 

Ahk = {uh\^ through the formula 



Tin (Ag(t)) =J2 Tr n (e* 

hk 



n Ht a ] h a k e-i Ht g) A hk 



(i) 



Tr HW (^'gU(t) 



(2.6) 
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One has to study the expression e^ m a h a k e~^ Ht , exploiting the fact that the expectation 
values (A) t are "slowly varying" if the matrix Ahk is "quasi-diagonal" (if Ahk = 5hk, A is 
a conserved charge). We give here only a sketchy account of the main points (for details 
see It proves useful to use in the Heisenberg picture a formalism in B(TL) reminiscent 
of usual scattering theory in TC, by means of superoperators, typically: 

= £[£,.], tc = z -[h q + h m ,.], V = jr\y,-]. 

In such a context operators of the form a h a k are "eigenstates" of TC with eigenvalues 
i (E h - E k ). Setting U'{t) = e w ' 4 one has: 



U\t)(a\a k ) = (lA\t)a{)(u\t)a k 



(2.7) 



+ioo+»j +ico+ri 
-ioo+?7 — ioo+r? 



1 



1 



7T + 



7rT(z) 



where T(-z) = V + V 



rV. 



(2i 



z — H' z — H z — H z — H z — Tt' 

T(z) is reminiscent of the usual T- matrix of scattering theory and plays a central role in 
this treatment: it will be called "scattering map". The operator T(z) has poles on the 
imaginary axis for z = \{e a — e@), e a being the eigenvalues of H. In the calculation of 
expression (|2.7Q we shall assume that the function T(z) for Rez w e with e ^> 5 (5 typical 
spacing between the poles) is smooth enough, so that the only relevant contribution from 
the singularities of (z — TC ) _1 stems from the singularities of (z — 7i ) -1 ; this smoothness 
property is linked to the fact that the set of poles of (z — Tt ) _1 goes over to a continuum 
if the confinement is removed yielding an analytic function with a cut along the imaginary 
axis, that can be continued across the cut without singularities if no absorption of the 
microsystem occurs. More precisely T(iy + e) is considered as practically constant for 
variations Ay w ^, where To has to be interpreted as a collision time. Treating expression 
( p.7| ) we make use of the inequality 



\Eh — E k \ -C — , 

T 



(2.9) 



whose physical meaning is that the typical variation time T\ of the quantities (A) t is 
much larger than r . One then arrives at the following very perspicuous structure for 
U'{t){a{a h ): 

U (t) {a{a k ) = a h a k + tC' {a h a k ) 



£ (a h a k ) 



% r 

h 



H eB ,a h a k - - ( f (1 \al 



a k — a h 



n 



E^S^S, (2.10) 



where H cS = Ho + V cti and 
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p(l) 

-"-fcA 



XX' Z 



XX' 

gr 



T { - % -E r + e\ a^j a) g + a r (t (^E g + e)ti g 



T ( -^£ r + e) aA a\ - a r [T (J-E g + e)a g 



2 ~ 3 v <Ai(r(-i^ + £ )a,)at|AQ 

2^ ^ , B B B 7^\ A \ a 9 



0* 



? E„ + E x - E k - E y -ihe 



\X')(X'\a g 



|A')<A'|a s 



and | A) denotes an eigenvector of _£/m with eigenvalue E\ and of H with eigenvalue zero. 
Eq. ( p.!0|) shows a typical structure arising in the calculation, which we will also find in 
the more complex situation examined in §3, where the form of the different operators is 
further commented on. Let us observe that V c r a g and T rg are not c-number coefficients, 
but operators acting in the Fock-space for the macrosystem, as stressed by the hats; they 
are connected respectively to the self-adjoint and anti-self-adjoint part of what can be 
considered as an operator valued T-matrix. The last contribution displays the "bilinear 



structure" of the third term in the r.h.s. of ( 1.2|) , connected to irreversibility and CP 
and not reproducible in the Hilbert space formalism, even resorting to an interaction 
potential which is not self-adjoint. Within the approximation leading to ( |2.10|) one has 
~ ^J2hx^hx^h\ an d therefore C' N = 0. Appealing to Q2.6Q we may obtain an 
evolution equation for the matrix elements Qf g which is meaningful on a time scale much 
longer than the correlation time for Sm'- 



dg, 



'gf 



dt 



jr (E g - Ef) Qgf + jrY, QghQhf + ^ E QgkQkf + J- E ( l ^)gk Qkh (L 



K) fh 



with 



Q 



kf 



H 



T ( ~j^E k + e J a k J a} f g M {t) 



T [--E k + e)a k )d) 



1 



E k + Ex — E 



f 



M 



ihe 



is a complete system of eigenvectors of Qu{t), (£m(£) = J2m) 7r C(t)l?(^))(^(^)D- To 



show the connection with (|1.2| ) we introduce in TiP^ the operators Q*' \ L^i : 



(*|Q l/> = Q 



kf , 



(k\txl\f) 



thus attaining in the Schrodinger picture the full evolution of g^\ given by the typical 
Lindblad generator: 



dgU 
dt 



Heff; Q 



(1) 



+ 



2h 



Q (1) + Q (1) t 



+ ± v[ (1) #)[ w t 



(1)4 



ft 
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where f\[J — H ^ + ~ ^Q*" ^ — Q*" Furthermore according to preservation of trace we 



A (i) A (i) f f (i) tf (i) 
have Q + Q T = - £ A? L a? T L Af ■ 



As it is well-known (|2.11| ) is apt to describe very different physical situations. If the 
last contribution, which we will call "incoherent", may be neglected, at least as a first 
approximation, eq. (|2.11|) is equivalent to a Schrodinger equation with a possibly complex 
potential. In the case of a particle interacting with matter this equation is well-suited 
to describe a coherent optical behaviour, for example in terms of a refractive index, as 
it is usually done in neutron optics ||, [TO]] and recently also in atom optics ]TT], [TJ]. In 



this framework the operator g/ ^ is to be interpreted as an optical potential, which in 
our formalism is naturally linked to matrix elements of the T-operator, thus showing the 
connection between the effective, macroscopic description through an index of refraction 
and quantities characterising the local interactions. The T-operator may be replaced 
by phenomenological expressions (for example the Fermi pseudo-potential in the case of 
neutron optics). This picture is particularly useful to deal with particle interferometry. 
To see how the last contribution may be linked to an interaction having a measuring 
character let us introduce the reversible mappings A tH , = U§1 ■ 17$, t, where U$ t , = 

Texp(— | J*, dt (Hq \t) + zQ \t))), corresponding to a coherent contractive evolution of 

the microsystem during the time interval [t',t"], and the CP mappings £ A £ = L A ^(£) ' 

L A ^t(t), having a measuring character, as it is clear from their very structure, reminiscent 

of the reduction postulate. The expression of the operators L A t shows how these mappings 
may be linked with a transition inside the macrosystem specified by the pair of indexes 
£, A, as a result of scattering with the microsystem. These transitions are in general not 
detectable, but under suitable conditions they could prime real events. The solution of 
( [2. 11| ) may be written in the form: 



Lh = AttoQto + [ dtiAtt 1 C\ 1 t 1 {ti)At 1 t Qt + • • • , (2.12) 



that is a sum over subcollections corresponding to the realization of no event, one event 
and so on. The set of variables N\^(t), t > to, (number of transitions up to time 
r), define a multicomponent classical stochastic process, and ( |2.12|) corresponds to the 



decomposition of the evolution map on the space of trajectories for N\^(t). This is a 
straightforward generalization of the typical "counting process" considered by Srinivas and 
Davies [13R. Of course other decompositions in terms of operation valued maps are possible 



on trajectory spaces related to different observables, as indicated in §1. When the first 
term in Q2.12j ) is largely predominant a wavelike description as given by the Schrodinger 
equation is sufficiently accurate and small disturbances, conveyed by the other terms, 
play no significant role. In a different physical context however, as would be the case for 
the brownian motion of a particle interacting with an ideal gas, the interplay between the 
contractive and the incoherent part has a major role. Being interested in the dynamics far 
away from the walls the quantum number h corresponds to the momentum variable p h , 
and supposing that momentum transfers are small (Fokker-Planck approximation) one 
arrives at an equation describing diffusion in phase-space which has the following form 
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sec 



0): 



dg 
~dt 



Heir, Q - D pp [x, [x, q]] - D qq [p, [p, 



17] 

2M 



% {p, q}] 



D qq , D qq and 77 being diffusion coefficients linked in different ways to the operators 
and M the mass of the particle. 



3 Theory of a macrosystem: thermodynamic evolu- 
tion by a scattering map 

We consider a very schematic model of macrosystem in the non-relativistic case built by 
one type of molecules with mass m confined inside a region u, interacting by a two body 
potential V(|x — y|); for the sake of simplicity no internal structure of the molecules is 
taken into account. In the field theoretical language the system is described by a quantum 
Schrodinger field (QSF): 



to 



£>/(x)a/ 
/ 



>fg 



(3.1) 



h 2 



-^ A 2«/(x) = E f Uf(x), w/(x) =0 x G du>. 



We shall assume the following Hamiltonian to take local interactions and confinement 
into account: 

h = J2 E f h Y a s + h V a h v hi2f2h a h" a h ( 3 - 2 ) 



hh 

V hhhh = [ d\f dV<(xX(y)y(|x-y|) M/2 (y) U/l (x). (3.3) 

J UJ J UJ 

Eq.( p.l|) is linked to the basic "local" Hamiltonian for the non-confined field (NC) 

n 2 



H NC = J grad^N C ( x ) " grad^ 



NC 



^ NC (x),^ c (x') 



V{r)ip 



5(x - x) 



N < I X + - ) ^ NC I X 



simply selecting the part of H NC related to the "normal modes" u f typical of the confine- 
ment: in fact the preparation procedure should imply a kind of relaxation of , 0nc( x ) to 
^(x). Skipping this problem and also the related question of the full explicit structure of 
£ for a realistic system, we shall simply take the Hamiltonian Q3.2|) containing only the 
normal modes of the field inside u. If we are interested in a hydrodynamic description 



8 



relevant observables are constructed starting with the densities of the typical constants 
of motion, mass and energy: 



m (x) = m f(x)^(x) (3.4) 

h 2 

e(x) = grad^fx) • grad^(x) 

2m 



\ L A * ( x - * ( x + D v(r) * ( x + D * ( x - D 



where the dependence of u x on x is generally negligible if V(r) is a short range potential. 
In the case of a kinetic description we replace (|3.4|) by the "Boltzmann" operator density 



/(x, p) = mJ2hk (x,p)|wfe)ofe, where F is the density of joint one particle 



position- momentum observables |T5|, [16|]. These densities lead to slowly varying quantities 
if they are integrated over regions large enough in space or phase-space, since one has 
constants if the integration is extended over the whole space. The constants of motion 
leading to this subdynamics are linked to very fundamental symmetries: time translation 
invariance and gauge symmetry. Our relevant observables have the general structure: 

^ra} h A hk (£)a k , ^ a} hi a\ 2 A hlh2k2kl (y.)a k2 a kl (3.5) 

hk k 1 k 2 



A hlh2k2kl (x) = - £ eft u* hi (x - u* h2 (x + V {r)u k2 (x + u kl (x - | 
We thus have to study in Heisenberg picture the expressions: 



^2en m a{a k e n m A hk {£,), £ en m a hl al 2 a k2 a kl e n m A hlh2k2kl (x.), (3.6) 

hk 



and shall take into account that by the slow variability only terms that are "diagonal 
enough" are really relevant; the sums should be restricted to indexes such that: 

— | E h — E k | < — — | E hl + E h2 — E kl — E k2 1 < — , 
n T\ n T\ 

where T\ is the characteristic variation time of the relevant quantities. We would like to 
stress the fact that the QSF is the basic tool to describe a massive continuum, just like 
the quantum electromagnetic field describes a massless continuum. The dynamics of the 
QSF, "0(x, t) = e^ Ht t(j(x)e~n Ht , in terms of which one can rewrite (3^3), is given by the 
simple field equation 

ih^i>(x,t) = ~A^(x,t) + J <fyft(y,t)V(\x-y\)$(y,t)#(x,t), (3.7) 

however no such equation holds for the expectation value of the field i>(x.,i) = (^(x, t)) 
due to correlations in the non-linear term; ^(x, t) is not useful to calculate the expectations 
of operators ( |3.6|) and therefore a classical Schrodinger field equation for ?/>(x, t) has no 
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physical meaning in general. In this respect the case of electromagnetism, where no self- 
interaction of the field occurs, is deeply different and allows classical electrodynamics 
to play an important role. To the macrosystem one associates typical "thermodynamic 
state" parameters: the velocity field v(x, t), the temperature field /3(x, t), the chemical 
potential field /u(x, t) in the case of the hydrodynamic description or more generally a 
field /i(x, p,t) on the one-particle phase-space in the kinetic case |T7[. The parameters 



/?(x, t) and //(x, t) (//(x, p,t)) determine the expectation values of energy density and 
mass density (the Boltzmann operator); let us briefly recall how the relation between 
state variables and expectation values is established [ 18(1 . At any time t one considers the 



whole set of statistical operators {w} which yield the expectation values assigned at that 
time: 

(e(°)(x)>, = Tr(e(°>(x)u,) , (p m (x)) t = Tr(p m (xH , (/W>(x, p)) t = Tr(/C )(x, p)w) 



g(°)( x ) = J- (ih^-mv(x,t))ft(x)- f-ift^-mv(x,t)Wx) (3.8) 



+ \ L A * ( x - ^ ( x + v{r) ^ ( x + * ( x - 1 

p^(x) = m^(x)^(x), / (0) (x, p) = /(x, p - mv(x, t)) 



where the quantities indexed by (0) (depending explicitly on the velocity field v(x, t)) 
represent densities in the reference frame in which the continuum is locally at rest. The 
velocity field is related to the expectation value of the momentum density p(x) through 
the relation (p^(x)) = 0, where 



p<°>(x) = ~[^(x) (-ih-^- - mv(x,t)) 



'. d \ % 

ih— mvfx, t) ij) + (x) 

. 9x / 



^(x) 



or equivalently (p(x)) i = v(x, t){p m (x)) t . Then one looks for a statistical operator in the 
set {w} such that the Von- Neumann entropy S = — A;Tr (w log to) is maximal, i.e. the 
most unbiased choice of a statistical operator leading to the given expectation values. 
The unique solution of this problem is 

e - J u d\ /3(x,t) [g(o) ( x )- At (x,t) / 5 m (x)] 

w[P(t)^(t)Mt)) = Tre _ L , M(x , t)[g( „ )(x) _ Mx ^ m(x)] (3-9) 
and analogously in the kinetic case. 

The corresponding S = — kTr (w[/3(t), p(t), v(t)] logw)[/3(t), fJ,(t), v(i)]) is the thermody- 
namic entropy of the macrosystem. If the time evolution of the expectation values 
(e^(x)) 4 , (p m (x))i, ((/(x, p))tj is given by the Hamiltonian evolution ( |3.6| ) or more gener- 
ally by a map M.' ttQ , having a preadjoint M.tt which does not decrease the Von-Neumann 
entropy, one immediately has that the thermodynamic entropy is non-decreasing. In this 
way one establishes the second principle of thermodynamics on a very clear dynamical 
basis. 
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In the simplest scheme of macroscopic dynamics the thermodynamic state parameters 
v(x, t), /3(x, t), /i(x, t) (/i(x, p, £)) at time to determine its evolution for t > to, e.g., by 
differential equations. Phenomenology shows that this is very often the case. Tackling 
the problem from the theoretical viewpoint one is induced, considering the operators 

e(x) = i[ff,e(x)] (/(x,p) = i[F,/(x,p)]), 



to calculate their expectations with the statistical operator given by (|3.9|) . This leads 
to wrong results as can be seen from the fact that the expectation values of the cur- 
rents which can be associated, through a conservation equation, to these operators would 



vanish |[19|| , due to time reversal invariance of microphysics, thus failing to describe any 
dissipative flow (e.g., heat conduction, viscosity, etc.). The idea of a time scale for the 
thermodynamic evolution and of a related subdynamics for the basic densities leads to 
a refinement of the aforementioned procedure: assume that %[H, ■] can be replaced by 
a mapping C' , initially defined on the linearly independent elements a} h a k , ala^a^a^, 
giving the slow time evolution of the relevant variables. In this way not only the sta- 
tistical operator w[/3(t), //(t), v(t)], but also the evolution operator is tuned to the rele- 
vant observables. Then one has the following set of closed evolution equations for the 
thermodynamic fields v(x, t), /?(x, t), ;u(x,t), (yu(x, p,t)) related to the basic observables 
A = p m (x), p(x), e(x), (/(x,p)): 

jTr(Aw[(3(t),f,(t)Mt)]) =Tr((£'A)w\J3(t),rit)Mt)])- (3-10) 

The non-Hamiltonian form of the map C' eliminates the aforementioned difficulties with 
vanishing dissipative flows; preliminary investigations of the consequences of (|3.10| ) in the 



case of a dilute gas indicate that it could be the right solution. The map C that ade- 
quately replaces t[H, ■] for the slow variables must generate an evolution of the relevant 

observables that preserves their positivity properties (e.g., p m (x), j(x, p)) and also con- 
servation of mass (M = J G&p^x) = J d 3 x<i 3 p /(x, p)) and of energy (E — j d 3 xe(x)). 
Then £ M = and C E = 0, while positivity with respect to observables constructed in 
terms of creation and annihilation operators could arise by a stronger property, reminding 
CP: 



)(ip h \U' (a} h a k ) ip k ) > (VV/J^' (ah^h^ak^) ^fe 1 fe 2 ) > 0, (3.11 



hk 



for any choice of {iph} and {iph^}- 
The time evolution of the typical expressions (|3.6|) can be studied by a procedure quite 



similar to that already shown in §2, based on the representation ( 2.7|) in terms of the 
"superoperator" Tio, having a} h cik, aLat a^a^ as eigenstates and of the superoperator 
T(z), which was called scattering map. If suitable smoothness properties of T(z) occur, 
essentially only the poles of (z — 7Y ) _1 contribute to the calculation of (|2.7| ), so that the 
following asymptotic representation holds: 

h 

U{t) (a{a k ) = a\a k + tC {a\a k ) r < t < — — ; (3.12) 

\hih — -fri- 



ll 



r is linked to smoothness properties of T(z) and can be interpreted as the typical duration 
of a collision between two particles interacting through the potential V(|x — y|); r fixes 
a time scale that is assumed to be much smaller than the typical variation time T\ of our 
relevant observables. £ is a linear mapping defined initially on the linearly independent 
elements a ] h ak, a^a^a^a^. For brevity we simply describe the structure of £ , skipping 
the derivation. The formalism produces the typical structure of two-particle QM, with 
an N-body correction due to the Pauli principle. A two-particle scattering operator is 
defined by 

f (2) (2) = y(2) + y(2) 1 yf tf W = £00 + V £\ ( 3 .13) 

where these operators, labelled by the index (2), are defined in the Hilbert space of 
two identical particles by matrix elements in the two-particle (symmetric or antisymmet- 
ric) basis [hh)] the matrix elements are: 

(Wil^l/a/i) = (E h +E f2 )±(5 hf2 5 hfl ±6 l2fl 6 hh ) 

(kklV^lhh) = V hl2hfl (3.14) 

(WilvfWi) = (l±h h ±n h )V hhhfl (3.15) 

the coefficients Ef, V^/^ are given in ( |3.2| ) and (|3.3[), the factor (l±n^ ±n; 2 ) is given in 
a more indirect way: the "two-particle" QM expressed by the aforementioned operators 
provides c-number coefficients in Fock-space operator expressions initially defined on the 
Fock-space basis \ . . .rif . . .), rif being the occupation numbers of the different field modes 
/: rif G N in the Bose case, rif = 0,1 in the Fermi case; therefore the factor (1 in^ ±n; 2 ) 
depends on the Fock-space basis elements on which the final Fock-space operator is acting. 
Also the adjoint operator in the two-particle Hilbert space will be useful: 



V {i) + V% (3.16) 



o-(2) 



The superoperator d consists of an Hamiltonian part ^[H eS , •] and of a part, analogous 
to the one in fl2.10Q , reminding the Lindblad structure ( |1.2| ). The formally self-adjoint 
Hamilton operator H cS is initially defined on the Fock-space basis \ ...rif...) by the 
following expression: 

H ctt = ^Efa}a f + -J2H 1 & i2 V hhf2fx & f2 d 'h ( 3 - 17 ) 



By comparison with ( |3.1| ) one can notice that introducing the time scale t ^> tq the 
coefficients Vi^f^ related to the basic interaction between the field modes is replaced by 
y\[i % M y i linked with a full, Pauli principle corrected description of the two body collisions 
in the medium, expressed in terms of the self-adjoint part of the operator f( 2 ) {£) . The anti- 
self-adjoint part ~ not zero if one goes beyond Born approximation 
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and provides a contribution to C! analogous to the second term in the l.h.s. of ( |iT2|) , 
of the form — | ( T^ 2 \a h — ajj T^ 2 \a k ), that due to sign "-" cannot be rewritten as 

f^ 2 ), • ; the operator f ( 2 ) is defined on the Fock-space basis by 



o E a h a h 

Z hf2 



* (l 2 l 1 \-{rV)(E h +E f2 +ihe) 



f^(E h + E h +ihe) ')\hh)a f2 a h 



At this point one immediately expects a third contribution to d related to the product 
structure a} h a k and involving both f( 2 ) and T^ 2 ) ; this contribution is given by h Y^x-R-hX^ -^fcA 
and reminds the structure of the third term at the l.h.s. of the Lindblad expression ( |1.2|) , 
where it accounted for decoherence (or state reduction, or event production). The oper- 

~ (2) 

ators R kX are defined on the Fock-space basis by: 



R 



(2) 
kX 



-i^2e(l±h x ±h k ) J2 



hh 



{k\\J^{E h + E h +ihe)\f 2 f 1 ) 
E k + E\ — E h - E h - ihe 



a h a h 



The factor J 2e(l ± h\ ± h k ) arises in the approximate factorisation of a Pauli correction 
term depending both on h k and hh'- 



2s ( 1 ± h x ± -{h h + h k 



2e(l ± h x ± h h )J2s(l ±h x ± h k ), 



(3.18) 



this factorisation, which is a good approximation if the Pauli corrections are not very 
large, is a typical quantum condition, which together with r <C T\ must be satisfied 
for the validity of the simple thermodynamic behaviour that we are considering in this 
section. The final structure of L is formally the same as in ( |2.10| ) 



C'a f h a k = % - \H eS , ala k ] - - ( [f (2) , 
n 1 J n. vl 



a k - a' h 



f(2) 



n 



E^Ml (3-19) 



the main difference lying in the space in which these operators act, according to the two 
different physical situations. As a consequence of unitarity of U one can prove that within 
the approximation leading to expression ( |3.19|) one has: 



f (2) 



4 



hX ) 



(3.20) 



hX 



therefore one can replace the expression f^ 2 ^ « iX^a^Ia^Ia * n dOH) : i n this way 
the conservation relation C' M = is exactly satisfied. It can be easily shown that 
Shfc('0/i|([l + r £ Wh^k)ipk) > to first order in r, so that the positivity property (|3.11|) 
is satisfied. 

As a preliminary check of the formalism let us report that the calculation of L (a h a k ) 
yields the typical structure of the collision term of the Boltzmann equation with the Pauli 
principle corrections: Jflgt^g and — i r^,a}j + c.c. are respectively the "gain" 

and the "loss" part of the collision term; C' (a k a k ) yields also the streaming term of the 
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Boltzmann equation. The study of C (cl^ a^a^) is not yet finished, it should end 
up with a full hydrodynamic and kinetic description of a one component continuum. 
The approximations leading to £ are based on the smoothness assumption related to 
the condition r <C T\ and to a "one mode" approximation for the description of the 
dynamics in the time interval r . A finite parameter e ~ Ht appears in the formal 
expression of C ; the final results do not appreciably depend on this by-product of the 
approximations if To -C T\. In a sense the simple thermodynamic behaviour expressed by 
( p.lOj ) arises by an approximation and this is indicated by the presence of e: an appreciable 
dependence of the results on e indicates a failure of the smoothness assumption and 
of the related approximations. Let us stress finally that in this approach existence of 
closed evolution equations for the thermodynamic state variables avoids any factorisation 
assumption for the distribution functions and therefore goes far beyond the approach 
based on the truncation of a hierarchy. 



4 Dynamics with memory effects 

According to § 3, when t[H, ■} can be replaced by the map C' given by ( 3.19| ) the family of 
generalised Gibbs states w(t) replaces the family of statistical operators g t — U(t — t )g to 
(for simplicity we assume now time-translation invariance) . This means that to determine 
the evolution of the thermodynamic state from a given time t onwards nothing else than 
(3, fx, v at that time has to be taken into account: i.e. no bias comes by the previous 
history f3(x, t'), /i(x, t'), v(x, t'), t' < t. This is no longer true if the condition t\ r [see 
eq. (|3.12p 1 is not satisfied. Let us assume that at an initial time T the statistical operator 
qt can be identified with the Gibbs state w(T) related to it: 

Qt = w(T). (4.21) 



By a straightforward calculation [18| one has: 

. . -</3(T).e(0)h(t-T)]) + ([MT)/3(T)].p m h(t-T)]> 

p t = e -W~ T ) o T e^ H{t - T) = — (4 221 

Ut UT Tre -{/3(T).e(0)[-(t-T)])+([ / ,(T)/3(T)].p m [-(t-T)]> ^ " > 



Tr e 



-</3(t)-e(°)>+<[ M (t)/3(t)].p m >- J* dr^(m-T)-m(-T))-Mt-T)P(t-T)]-p m (-r))) 



where A(r) = e^ HT Ae~^ HT and (f3(t) ■ A) = J^d 3 x/5(t, x)A(x). In the last term of the 
exponent the history of the thermodynamic state during the time interval [T, t] appears; 

by Pm = — divJ m , e = — divJ; it can be rewritten in the more perspicuous form: 

i_T dr L (</?(* - r) ■ e(°)(-r)) - ([»(t - r)f3(t - r)] • p m (-r))) = (4.23) 



dr 

t 



dr' 



T 



(^(3(r f ) ■ e<V " 0> " (^(r'Hr')} • fflV - t)) 



- (grad^(r') • jj 0) (r' " t)) + (grad(/3(r')^(r')) • J^(r' - t)) 
+ £ dr' da n • [{(3{r\ x) J® (x, r'-t))- (/3(r', x)/i(r', x) (x, r' - t)) 



14 



where time and space derivatives of the thermodynamic fields appear on the same footing; 
by the last term also matter and energy exchanges with the environment during the 
preparation time [T, t ] can be described. Taking expression ( |4.22| ) for g t to calculate 
the basic expectation values and thus determining the thermodynamic state variables 
for t > t as in (|3.10|) , one has for them closed evolution equations having as input 
/3(x, r'), //(x, r'), v(x, r') r' G [T, to). Such equations are generally used and work under 
the hypothesis that memory decays within a typical correlation time. This also helps 
to attenuate the problem of the initial choice Q4.21]) . We mention in passing that no 
problem about condition ( 4.21|) exists if one assumes the point of view of "informational 
thermodynamics": then qt = w(T) is just dictated by the measured values of the relevant 
variables at time T; however this approach does not explain why the previous history of 
a concrete collection of macrosystems is irrelevant just before T. Let us also mention 
the solution given by Zubarev: T is shifted to — oo, thus taking off any previous history; 
however this limit is highly critical and since also a thermodynamic limit is involved, it 
shifts the problem of thermodynamic evolution to a cosmological one. In our framework a 
time scale is associated to the system and the generator £ should have a non-Hamiltonian 
part. This provides a mechanism by which memory can decay. Let us assume that the 
preparation time t — T is larger than the decay time of the memory: at this point the 
history that comes before T is irrelevant for the dynamics that C is able to describe. Then 
the choice ( |4.21| ), that is not biased by this history, is adequate. When C = l[H, ■] + C , 



calling q([/3, fi, v], t, T) the operator on the l.h.s. of ( |4.22| ) one has 



Tr (Au(t, T)qt) = Tr [(W>^ (r) i) q(\J3, fi, v], t, T) 



n 
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